Nuclear binding, correlations and the origin of EMC effect 
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Recent data for the slope of the EMC-ratio in the intermediate x-region for light nuclei, with 
3 < A < 12, have the potential to shed new light on the origin of the EMC effect. Here we study 
the role of nuclear binding using the scaling variable y, best suited to take into account this effect, 
and the understanding of the average nucleon removal energies, E, provided by state-of-the-art 
calculations based on nuclear many body theory. We find an excellent correlation between the new 
EMC data at x ~ 0.5 and E for nuclei with A from 3 to oo, indicating that in this x region binding is 
an important ingredient to explain the EMC effect. The role played by nucleon-nucleon correlations 
in this context is also discussed. 

PACS numbers: 25.30.Fj, 24.85+p, 13.60Hb, 21.60De 
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The so called EMC effect was, in 1983, a surprising 
discovery [![. The ratio R(x) between the Deep Inelastic 
Scattering (DIS) cross sections of leptons on iron and the 
deuteron was found to sizably differ from unity. At low 
values of the Bjorken scaling variable x the ratio is < 1, 
at medium values of x R drops from ~ 1 to values as 
low as 0.8 and at large x it reaches values larger than 1. 
While the latter feature can be quantitatively explained 
by the smearing of the parton distribution functions aris- 
ing from the momentum distribution of nucleons in nu- 
clei, the former one is accounted for by including the 
effect of nuclear shadowing. The reduction of R(x) at 
intermediate x was not easily explained, although many 
different models have been proposed (for reviews of the 
extensive literature see, e.g., Refs. 0-3 )■ Often, but not 
always, these models do involve the role of the binding 
of nucleons in the nucleus. In general, however, the ef- 
fect of binding alone is not large enough to reproduce the 
data In particular, the rise towards small x seems to 
require the contribution of pions — enhanced due to nu- 
cleon binding — Q and consideration of the role of the 
recoil final state interaction. No "definitive" explanation 
has surfaced, in part due to the fact that no model pro- 
vides a quantitative description of the whole range of x, 
unless very specific assumptions are made. 

Recently, new measurements have been performed by 
Seely et al. Q with the hope to elucidate things better by 
adding precise data on light nuclei. For these nuclei, the 
nuclear structure is extremely well understood, mainly 
thanks to Green's Function Monte Carlo (GFMC) cal- 
culations based on state-of-the-art models of the nuclear 
hamiltonian, strongly constrained by the observed prop- 
erties of the two- and three-nucleon systems. For light 
nuclei, nuclear properties change rapidly as a function 
of the mass number A, so that different models for the 



EMC effect could yield rather different predictions. The 
data of Ref. 0] mainly cover the intermediate region of 
0.35 < x < 0.7, which is little affected by smearing and 
away from the low- a: region, where the potential contri- 
butions of pions, coherence effects and final state inter- 
actions complicate matters. 

In order to provide data least affected by experimental 
errors, the authors of Ref. 0] have studied in particular 
the slope dR/dx, which is insensitive to normalization er- 
rors. The correlation between dR/dx for A=3,4,9,12 and 
the average density has been studied in order to find 
a phenomenological relation that could further elucidate 
the EMC-cffect at mid- a;; no unambiguous (ilinear) cor- 
relation has been found. In particular, the dR/dx value 
for 9 Be does not follow the tendency observed for the 
other light nuclei. The relation to quantities such as nu- 
clear binding energy or nucleon separation energy Q also 
has produced no further insight. A correlation with the 
inclusive quasi-elastic cross section for electron scattering 
at x > 1 has been found 0, [l(| , and attributed to the 
role of nuclear high-momentum components that increase 
with increasing A. However, in the standard treatment 
of Fermi motion, folding the nucleon parton distributions 
with nuclear momentum distributions alone would lead 
to values of the EMC ratios > 1 and a mid- a: slope dR/dx 
of the wrong sign. 

In this paper, we argue that DIS data are best ana- 
lyzed in terms of the scaling variable y, widely employed 
in studies of a variety of scattering processes involving 
composite targets. We explore the correlation of dR/dy 
with nuclear binding, which in many approaches is ad- 
vocated as an important element for the explanation of 
the EMC effect, and is responsible in particular for the 
decrease of the EMC-ratios to values < 1 at intermediate 
values of the scaling variable. 
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In general, DIS cross sections are studied in the frame- 
work of the parton model, which explains the occur- 
rence of scaling and provides the basis for the deter- 
mination of the parton distribution functions 11]. The 



constituents the lepton scatters from are assumed to be 
point-like spin one half particles carrying, in the Infi- 
nite Momentum Frame (IMF), a fraction x of the nu- 
cleon four-momentum. This quantity is identified with 
the Bjorken scaling variable x — Q 2 / (2 Mi/), where Q 2 , 
M and v are the squared four-momentum transfer, the 
nucleon mass and the lepton energy loss, respectively. In 
the IMF interaction effects are expected to vanish [l2j ]. 
and the partons in the nucleon are on-shell. Within this 
approach the binding (off-shellness) of nucleons in nuclei 
is conceptually difficult to introduce. 

In order to avoid this problem, we here employ a some- 
what different approach, first developed to analyze quasi- 
elastic electron-nucleus scattering, a process perfectly 
analogous to DIS. Scaling of the quasi-elastic scattering 
cross sections has been studied in terms of the scaling 
variable y [3] which is derived from the kinematics of 
the underlying process of elastic scattering from nucle- 
ons bound in nuclei. 

Energy-momentum conservation of elastic scattering 
from a constituent of mass to k yields 



(1) 



with q being the 3-momentum transfer, and to being, 
for quasi-elastic scattering, the mass of the nucleon, i.e. 
rn = M. Note that in this approach the initially bound 
nucleon is assumed to be off-shell. Off-shell means that 
there is no — \/|k| 2 + to 2 term on the right-hand side, 
only a term to; while the nucleon has momentum k it 
is assumed to be bound with total energy =0. As usual, 
the interaction energy of the knocked-out nucleon with 
momentum k + q is neglected, which at large |k + q| in 
general is a good approximation. An additional term E > 
0, accounting for the average nucleon removal energy, is 
normally added on the right-hand side of Eq.([T]). 

Neglecting the component of k perpendicular to q, 
which is justified in the limit of large |q| = q, yields 



the standard scaling variable introduced in Ref. [13 



y 



\J v 1 + 2mv — q . 



(2) 



The physical meaning of y is straightforward: it is the 
component of the momentum of the initially bound nu- 
cleon parallel to q in the rest frame of the nucleus. 

Here, we specialize y to the conditions appropriate 
for DIS on the nucleon, i.e. scattering from (basi- 
cally) up- and down-quarks, the rest masses of which, 
to = M u ,Md <10 MeV, can be safely neglected at the 
energies relevant to DIS. In this case, the expression for 
y of Eq. ([2]) simplifies to 



where a 'tilde' has been added to remind ourselves that 
y corresponds to the to — )• limit 14 1. The physical 
meaning of y is analogous to the one of y in quasi-elastic 
electron-nucleus scattering: y is the component of the 
u/d-quark momentum parallel to q in the rest frame of 
the nucleon. 

Due to the large energy and momentum transfer in DIS 
on nuclei, the system of hit quark plus remaining debris 
of an initially bound nucleon will leave the nucleus with 
high energy and momentum, and the interaction with the 
(A-l)-nucleus is weak. The removal of this system costs 
an energy equal to the nucleon mean removal energy E 
which the scattered lepton has to provide. For the DIS 
process on the nucleon proper only the energy v' = v — E 
is available. This leads to 



y = V — q = v — E 



(4) 



As both the energy of the quark in the nucleon and E 
are well defined quantities in the rest frame of the nu- 
cleon (nucleus) they can be added without conceptual 
difficulties. 

For the isolated nucleon, i.e. for E = 0, using the 
variable y is known to leads to a scaling of the DIS cross 
section which is even better in quality than the scaling 
observed in terms of Bjorken x. This can be shown eas- 
ily, as y is trivially related to the Nachtmann variable 
£ [15[. While being usually defined using the much less 
transparent equation 



C = 2x/{l + ^l + AM 2 x 2 /Q 2 ) 



(5) 



y = v 



(3) 



Nachtmann's £ can easily be shown to reduce to (q — 
i/)/M, implying y = — £M. Using Nachtmann's £, which 
becomes identical to x in the Q 2 jv 2 — > limit, is known 
to extend the scaling property of DIS to lower v. 

As compared to x or £, y has a well defined physical 
meaning [36( in the nucleon rest frame, the coordinate 
system where the DIS experiments are usually done and 
where theoretical studies, such as lattice QCD calcula- 
tions, hope to provide a quantitative understanding of 
the nucleon structure functions 

In order to calculate E directly one would need to know 
the spectral function S(k,E), giving the probability to 
find in the nucleus a nucleon with momentum k and re- 
moval energy E. 

Theoretical calculations of the spectral function re- 
quire the knowledge of both the A-nucleon ground state 
and the full spectrum of eigenstates of the (A— l)-nucleon 
system. In addition, they involve an intrinsic degree of 
complexity, associated with the evaluation of non diago- 
nal nuclear matrix elements, that rapidly increases with 
A. 

Studies based on nuclear many body theory and real- 
istic hamiltonians have been carried out for A=3 T^- 1|| 
and for isospin symmetric nuclear matter [l9l [20j , while 
in the case of nuclei with A > 12 approximate spectral 
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functions obtained using the local density approximation 
are available [21| . In the present work, aimed at studying 
the correlation between binding effects and the EMC ef- 
fect for 3 < A < 12, we exploit Koltun's sum rule [22H23]] 
to obtain the average removal energies from the results 
of GFMC calculations. 

In the absence of three-nucleon interactions, Koltun's 
sum rule states that 

- (0) 



E 1 , - A 
— = - \T — 
A 2 A 



1 



E] 



where Eq/A is the nuclear binding energy per particle 
obtained from nuclear masses and includes a (small) cor- 
rection for the Coulomb energy, 



T = 



k 2 

d 3 kdE —S(k,E) 
2M v ' ' 



and 



E = / d 3 kdE E S(k,E) . 



(7) 



(8) 



The small contribution of the three-nucleon potential V3, 
which is known to be needed to achieve a precise deter- 



mination of Eq for A > 3, can be taken into account [24 1 
by adding a term (V^/A, where (...) denotes the ground 
state expectation value, to the right hand side of Eq.©. 

It has to be pointed out that the Koltun sum rule is 
an exact result, although its experimental verification in- 
volves severe difficulties. The results of the analysis of 
the (e, e'p) reaction carried out by Bernheim et al. 
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suggesting that the sum rule is badly violated, are due 
to the limited kinematical range covered by the experi- 
ment, which does not include the contribution of correla- 
tion effects, leading to the appearance of tails of S(k, E) 
extending to large energy and momentum (25j . 

The GFMC approach allows one to obtain essentially 
exact binding energies for the ground states and very 
good estimates of the energies of low-lying excited states. 
The wave functions resulting from GFMC calculations 
have also been employed to obtain density and momen- 
tum distributions, electromagnetic form factors and spec- 
troscopic factors, as well as to compute many electroweak 
processes of astrophysical interest. These studies have 
clearly shown that the ab initio approach based on the 
numerical solution of the many body Scrddinger equation 
and modern nucleon-nucleon interactions, fitted to NN 
scattering data, is capable to describe the full complexity 
of nuclear structure, including single particle properties, 
correlations and clustering. 

For all the light nuclei A < 12 of interest here, GFMC 
calculations have been carried out [for the lighter nuclei, 
exact solutions of the Schrodinger equation are also avail- 
able from simpler approaches, such as Variational Monte 
Carlo (VMC)]. From the resulting binding energies, mo- 
mentum distributions and (small) 3-body contributions 
we have calculated the average removal energies E which 
are listed in the following table. 



nucleus 


removal 
energy 


method 


reference 


2 H 


2.2 MeV 




[26] 


3 He 


14.6 MeV 


GFMC 


[27] 


4 He 


35.8 MeV 


VMC 


[27] 


9 Be 


43.8 MeV 


GFMC 


[28] 


12 C 


52.2 MeV 


GFMC 


[28] 


NM 


70.5 MeV 


FHNC 


[29] 



TABLE I: Average removal energies used in this work. 



The values of E are significantly larger than the typi- 
cal values that have been used in previous studies of the 
EMC effect [371, studies which often found that binding 
effects alone were not large enough to explain the EMC 
data. For 12 C, for example, the average removal energy 
one would derive from the centroids of the s- and p-shell 
removal energies as measured by (e, e'p) experiments [I(| 
is of the order 25 MeV. A similar value is obtained from 
mean-field calculations, such as Hartree-Fock. The dif- 
ference to the value listed in Table 1, 52 MeV, is due 
to the fact that the s- and p-shell peaks [13] account 
only for (the ~75% of) nucleons occupying mean-field 



orbitals [31] . The short-range correlations between nu- 
cleons, induced by the strong short-range components of 
the NN interaction (both central and tensor) lead to the 
appearance of nucleons in states of high momentum and 
high removal energy; the corresponding strength, which 
is thinly spread over a large range of initial momenta k 
and removal energies E, gives an important contribution 
to E. This picture has been confirmed by (e, e'p) experi- 
ments designed to provide a measurement of the spectral 
function S(k,E) at large values of k and E [25l |. 

Coming back to the EMC-effect, we show in Fig. 1 
the correlation between the average removal energy E 
and the derivatives dR/dy determined from the dR/dx 
values of Ref. Q [H|]. As Fig. 1 shows, there is an 
excellent linear correlation between E and dR/dy. This 
correlation is much better than the one of dR/dy (dR/dx) 
with other quantities (referred to above). This indicates 
that the driving quantity for the EMC-effect at mid- a; 
values is indeed the binding of nucleons, as described by 
the mean removal energy. 

To the extent that quasi-elastic electron-nucleus scat- 
tering at large x could be identified with high-momentum 
components, the correlation with ai found in [^, [l(| ac- 
tually would be based on the same physics as discussed 
here. There is a strong correlation between E and T (for 
12 C, for instance, £=25.2MeV, f=31.2MeV, the differ- 
ence leading to a comparatively small Eo/A=6AMeV). 
But while high-fc nucleons alone lead to EMC-ratios 
R > 1 and a positive slope near x ~ 0.5, the binding 
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E (MeV) 



FIG. 1: Correlation between the average nucleon removal en- 
ergy E and the slope of the EMC-ratio dR/dy obtained from 
the dR/dx values of 0,[32| at x ~ 0.5. Data points are labeled 
by the nuclear mass number A. 



leads to a much larger effect in the opposite direction, 
producing overall the R < 1 and dR/dx < 1 as found by 
experiment. 

In Figure 1 we have also included the point correspond- 
ing to uniform, isospin symmetric, nuclear matter (NM). 
The average removal energy has been determined as dis- 
cussed above from the variational results of Ref. [iij ]. 
obtained using the Fermi Hyper-Netted Chain (FHNC) 
summation scheme. The slope dR/dy has been extracted 



from Ref. [13] 0. 



As the determination of R{x) for nuclear matter in- 
volved a fit to the world EMC-data for all nuclei with 
mass number A > 12, the NM point is indicative of the 
fact that the excellent correlation between dR/dy and E 
is valid for all nuclei. 

We note that the correlation between dR/dx and E 
(not shown) is very similar to the one observed in Fig. 
1. Only the NM data point, which on average corre- 
sponds to larger Q 2 than the data of Q , would be slightly 
shifted due to a different conversion factor between x and 
y. While numerically the difference between x and y is 
small, these quantities differ radically in their physical 
interpretation. 

In principle, the approach based on many body theory 
employed in our work, while not including some of the 
mechanisms which are believed to determine the low- a; 
behavior of EMC ratio, may be used to obtain theoretical 
estimates of its slope at mid x. However, to achieve the 
level of accuracy required for a meaningful comparison 
with the data, one would need spectral functions com- 
puted using the Green's function Monte Carlo technique 
and including the full set of eigenstates of the recoiling 
nucleus, which are not yet available. 

In the case of infinite nuclear matter, S(k, E) has been 
computed within the FHNC/SOC summation scheme, 



including the contributions of one hole and two hole-one 
particle states However, compared to the calcula- 

tion of the ground state expectation value of the hamil- 
tonian discussed in Ref. [23], the work of Ref. [l9| was 
based on a oversimplified treatment of the three-nucleon 
interactions and involved a number of additional techni- 
cal difficulties (e.g. the orthogonalization of correlated 
states), leading to a larger theoretical uncertainty. As a 
result, the value of the average removal energy obtained 
from the spectral function of Ref. E « 61 MeV, 

appreciably differs from the one reported in Table 1. 

In conclusion, we have shown that there is a strong 
correlation between the EMC-effect at mid-x and the av- 
erage nucleon removal energy. This correlation covers all 
nuclei, from 3 He all to way to infinite nuclear matter. 
This confirms that the binding of the nucleons in the 
nucleus is very important for the EMC effect at x ~ 0.5. 

As binding plays an important role, our study of the 
data was done in terms of the scaling variable y. Be- 
ing derived as a property of initially off-shell quarks in 
the Lab frame, y can be generalized to take into account 
the additional off-shellness due to nuclear binding with- 
out conceptual difficulties. It should also be emphasized 
that, while y is particularly suited to discuss binding ef- 
fects, it may also be preferable to Bjorken-x in general as 
a scaling variable. Not only does y yield better scaling, 
it also has a more intuitive physics interpretation as the 
momentum component of the u,<i-quarks parallel to q in 
the rest frame of the nucleon, and allows for a unified 
description of inclusive scattering. The analogous occur- 
rence of y-scaling, and its straightforward interpretation, 
have been exploited to extract valuable information from 
the analysis of a variety of scattering processes, ranging 
from photon scattering from electrons bound in atoms 
to neutron scattering from quantum liquids and quasi- 
elastic electron- nucleus scattering (33[. 

We thank Robert B. Wiringa for providing the 
VMC/GFMC momentum distributions and three-body 
force contributions employed to obtain the results of Ta- 
ble 1. 
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